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Abstract 

The Casimir force due to a scalar field on a piston in a cylinder of radius r with a 
spherical cap of radius i? > r is computed numerically in the world-line approach. A 
geometrical subtraction scheme gives the finite interaction energy that determines the 
Casimir force. The spectral function of convex domains is obtained from a probability 
measure on convex surfaces that is induced by the Wiener measure on Brownian bridges 
the convex surfaces are the hulls of. The vacuum force on the piston by a scalar 
field satisfying Dirichlet boundary conditions is attractive in these geometries, but the 
strength and short-distance behavior of the force depends crucially on the shape of 
the piston casing. For a cylindrical casing with a hemispherical head, the force for 
a/i? ~ does not depend on the dimension of the casing and numerically approaches 
~ — 0.00326(4)/ic/a^. Semiclassically this asymptotic force is due to short, closed 
and non-periodic trajectories that reflect once off the piston near its periphery. The 
semiclassical estimate —he/ (9Q-Ka?)(\ + 2^ — r'^/a) for the force when a/r <^ r/R < 
1 reproduces the numerical results within statistical errors. 

PACS: 
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1 Introduction 



Until fairly recently, only the Casimir force between parallel plates ^4 |2j| could be reliably 
computed. Although the multiple scattering expansion of Balian, Bloch and Duplantier|3l 
H] allows for more complicated geometries, it was not used extensively to investigate the 
geometry dependence of Casimir forces beyond lowest order. Renormalization issues further 
obscured the relationship between finite Casimir energies and physical effects. Boyer[5], [6] 
calculated the (finite and apparently repulsive) Casimir pressure on an idealized metallic 
spherical shell in 1968 and DeRaad et al.|[7] the pressure on a cylindrical one in 1981. It 
later became clear that the finiteness of such Casimir self-energies depends rather sensitively 
on the assumed ultrathin ideal metallic boundaries [8], [9] and the physical interpretation of 
these results is debated. 

In the past few years this situation has changed dramatically. A better understanding of 
the relation between scattering theory and finite Casimir energies pjOj allows one to in princi- 
ple compute Casimir forces between disjoint bodies with arbitrary optical characteristics. A 
multipole expansion of the scattering matrix is particularly suited when the distance between 
the bodies is large whereas a semiclassical approximation becomes appropriate and rather 
accurate at short distances. The interaction vacuum energy of disjoint bodies is finite and 
the problem of determining the Casimir force between them in this sense has been largely 
solved - ignoring some issues at finite temperature [TT]. 

Below I define finite subtracted vacuum (interaction) energies that determine Casimir 
forces between bodies that are not separate. A similar geometrical approach was first used 
by Power [12] to obtain the well-known Casimir force between two parallel plates without 
regularizing potentially infinite zero-point energies. Svaiter and Svaiterp^ emphasized that 
these Casimir pistons show that certain (finite) physical vacuum effects do not depend on 
the regularization and renormalization procedure. Recently there has been renewed interest 
in such geometries [H] because the electromagnetic vacuum self-energy of a cube (as that 
of a sphere [5]) is positive [15]. The force on a partition in a parallelepiped nevertheless is 
attractive [in] at any position. For scalar fields satisfying Dirichlet boundary conditions this 
is a consequence of reflection positivityflUl fT7] . 

Semi-classical considerations suggest that the Casimir force on a piston strongly depends 
on the shape of the casing. The vacuum force due to a massless scalar field satisfying Dirichlet 
boundary conditions is here obtained numerically for some generalized Casimir pistons. The 
world-line approach to Casimir interaction energies of Gies et al. [18] is thereby generalized 
to geometries with connected boundaries. In Casimir pistons all domains are bounded and 
the mathematical treatment in fact is much simpler and clear-cut than for the scattering 
situation with a continuous spectrum. A finite subtracted vacuum energy gives rise to the 
force on the piston. In the examples studied here, all domains are convex and numerical 
computations are vastly simplified by considering the convex hulls of Brownian bridges. 
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2 World-line approach to the interaction energy of con- 
nected bounded domains 



Consider the heat kernel operator = e^^''^ for the Laplacian A with Dirichlet bound- 

ary conditions on a bounded domain C M'^. The spectrum of eigenvalues {A^ > 0,n G N} 
of the negative Laplace operator in this case is discrete, real and positive. The corresponding 
spectral function (or partition function or trace of the heat kernel), 

0^,(/?) = Tri^^,(/3) = 5^e-^^"/^ (1) 

nSN 

is finite for f3 > 0. It has the well-known [19], EHl EH [22] high-temperature (short-time) 
expansion, 

1 °° 

MP ~ 0) ~ J2(27rPr/'a4^) + 0{e-''l\ (2) 

For smoothly bounded domains, the Hadamard-Minakshisundaram-DeWitt-Seeley coeffi- 
cients a„(S)) in this series are integrals over powers of the local curvature and reflect av- 
erage geometric properties of the domain and its boundary [23l [2^ . For a bounded three- 
dimensional flat Euclidean domain 2), ao(S)) gives its volume and ai(2)) = — iSx)/4 gives 
the surface area of its boundary [23]. 02(2)) is proportional to the integrated curvature 
[sharp edges of the boundary also contribute [25j,i4j] and the dimensionless coefficient 03(2)) 
depends on topological characteristics of the domain [such as the connectivity of its boundary 
and the number and opening angles of its corners [23] H]]. Since 04(2)) 7^ implies a log- 
arithmic divergent vacuum energy that prevents a unique definition of the Casimir energy, 
this coefficient is crucial. Its geometric interpretation [23] is, however, not easily explained. 
Non-analytic and (for /5 ~ 0) exponentially suppressed contributions to the asymptotic ex- 
pansion of 0d(/3) are associated with classical periodic- and diffractive- orbits[2n| of minimal 
length /. 

The world-line approach to Casimir energies [18] is based on the fact[23l [27| that the 
spectral function for a bounded flat Euclidean domain S) can be expressed in terms of its 
support of standard Brownian bridges. In three dimensions, 

^»'« = X(2|fI^IWcS]. (3) 

where ^^(x) = {B^(x, /?),0 < r < /3;Bo(x, /5) = B^(x, /5) = x} is a standard Brownian 
bridge from x to x in "proper time" /5 and P[£^(x) C X)] denotes the probability that such 
a bridge is entirely within the bounded domain D. Note that V\tp{y^ C S}]/(27r/3)^/^ is 
the Green function from x to x in time /? of the associated diffusion problem with Dirichlet 
boundary conditions on 52). 

Eq.(l3]) implies that the spectral function for a domain of flnite volume is flnite. Di- 
vergences arise only in the corresponding zero-point energy. For flnite vacuum energies, 
the leading (flve, in three dimensions,) coefficients in the asymptotic power series of Eq.([2]) 
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would have to vanish. Ahhough impossible for a single domain, one can improve the asymp- 
totic behavior by considering a (finite) linear combination of spectral functions for domains 
{Dfc;A; = 0,l,...,M} 

0(/3) = ^c,0s,(/5) . (4) 

k 

When the coefficients Ck are such that 

^Cfca,(Dfc)=0 for i = 0,...,4, (5) 

k 

the "interaction" vacuum energy, 

is finite because the integrand is for /? ~ when Eq.dS]) holds. The £^aci!^k) ~ 

i Xln in Eq.Q are the (divergent) formal zero-point energies of a massless scalar 
field satisfying Dirichlet boundary conditions on the individual bounded domains T)k ■ The 
linear combination £^int of these vacuum energies may be interpreted as the difference in 
zero-point energy for domains of the same total volume, total surface area, local curvature, 
topology, etc... The subtractions have to be chosen so as to at most affect the physical 
quantity of interest in a calculable way. They thus could depend on the effect one wishes to 
describe. 

In the context of Casimir energies, a similar subtraction scheme was first used by Power [12] 
to calculate the original Casimir force between parallel metallic plates [Ij without regulators. 
Svaiter & Svaiter[T3] emphasized that the subtracted vacuum energy of parallelepiped pistons 
is physical and does not depend on calculational procedures. In reference [28] non-analytic 
contributions to the finite interaction vacuum energy were computed in leading semiclas- 
sical approximation. For numerical calculations, the subtraction method outlined above is 
preferable to conventional regularization procedures, because we shall see that it can be 
implemented by a restriction on the set of paths. One thus avoids the numerically difficult 
computation of small differences in large quantities. 

It is interesting (and will be useful numerically) to note that for convex domains 3D a 
Brownian bridge £/3(x) is wholly within S) only if its convex hull ^[£/3(x)] is entirely contained 
in S), that is 

V[lp{^) C D] = P[i3[£^(x)] C for S) convex . (7) 

In the continuum limit, a standard Brownian bridge process satisfies the stochastic dif- 
ferential equation [271 [29] . 

X — B fx B') 

rfB,(x, /5) = -^^^dr + dWr, with Bo(x, p) = x (8) 

(3-T 

where W,- is a (3-dimensional) standard Wiener process, i.e. each component of is 
normally distributed about the origin with variance r. Translational invariance and the 
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scaling property of random walks (and their convex hulls) imply that, 



VM^) C = P[£i(O)]0[x+ v^£i(0) C D] = P[£i(O)]0[x+ v^io[£i(0)] C . (9) 



Here multiplication and addition operations distribute over all members of a point-set, that 
is they scale and translate the whole loop or hull. Using Eq. to make the dependence on 
X and P explicit in Eq.(l3]) and inserting the result in Eq.(l6]) one obtains. 



int 



"8^ / V^^t^iW] CS),]),,(0) 



(10) 



where the expectation (. . .)e-^{o) is with respect to unit loops (standard Brownian bridges over 
unit time). Eq.l fTOi) is the basic formula used for computing Casimir interaction energies in 
the world-line approach[18]. Note that the overall proportionality constant in Eq. (fTO|) is 4 
times that in[18j because the unit loops here are defined by a standard Brownian bridge pro- 
cess rather than one of twice the variance [this can be verified by letting /3 — 2T in Eq. (fTO!) 
and comparing with[T8]]. Although the straightforward interpretation of an integrated in- 
teraction energy density is lost by the change in integration variables x = q/A,/3 = A"^, the 
last expression for Si^it in Eq. ffTOj) is slightly better adapted to numerical evaluation. 
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Figure 1: The interaction Casimir energy i?int(a) for a cylindrical cavity of radius r with a cap of 
radius R > r and a moveable planar piston at height a. i?int(o) is the difference in vacuum energy 
of the cavity with and without piston compared to this difference for a cylinder of the same radius 
r. This difference of differences in vacuum energies is finite for all values of a,r,L > 2r and R> r. 
i^int does not require regularization and for L ^ a only S^^-j depends on the position of the piston. 
Solid lines denote surfaces on which Dirichlet boundary conditions are imposed. Thick dashed lines 
in Figs, la) and I7) indicate where the surface of the cylinder would be: Brownian bridges that 
pierce these surfaces give no contribution to <Sint(a), leading to condition (fTTl) for large L ^ 2r. 
Note that all five different bounded domains are convex and bounded for finite L. 
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3 Algorithm for the Casimir pistons of Fig. 1 



Fig. 1 depicts a combination of convex domains that is suitable for computing the interaction 
vacuum energy Smt of a cyhndrical Casimir piston of radius r with a cap of radius R > r. 
A finite overall length L ensures that the volumes of all convex domains is finite. The 
restriction is of little concern here because the limit L — oo of Si^t exists and in fact is easier 
to calculate. Note that we could equally well have subtracted the vacuum energy at a fixed 
position of the piston (as in [121 1131 for the parallelepiped-piston), but taking the limit 
L — > oo then would have been numerically challenging. 

The main advantage of using the convex hull in numerical computations rather than the 
unit-loops, is that the hull is a vastly reduced point-set that retains all the information on the 
loop that is relevant for the partition function on convex domains [the sequential ordering 
and a great number of interior points of the loop being discarded]. The number of vertices of 
the convex hull increases only logarithmically with the number of points defining the loop. 
The convex hull of a loop given by a thousand points on average has ~ 55 vertices and that 
of one with a million points has about 190 vertices (see Fig. 2). The time required by efficient 
algorithms to compute the hull (with v vertices) of n points is proportional to nlog(f ) [31j. 
In evaluating the integrations in Eq. lfTOj) for sufficiently complex convex domains D like 
those in Fig. 1, one only requires the convex hulls. For high precision, i.e. for loops with a 
large number of points, n, it thus is advantageous to store the hulls rather than the loops 
themselves. It would be preferable to generate convex surfaces directly with the appropriate 
measure, rather than to construct them as hulls of Brownian bridges. But lacking a good 
algorithm to generate convex surfaces with the correct measure, they were here constructed 
as hulls of Brownian bridges [systematic errors due to the finite number of points of the 
original loops severely limit the accuracy one can achieve with this procedure]. 

To evaluate the spatial- and scale- integrals in Eq. (fTO!) for a given point-set i3[^i(0)] = 
{vj = (xj. Hi, Zi);i = 1, . . . , v}, note that (for r) the only significant contribution to the 
interaction energy £^int of the combination of domains shown in Fig. 1 comes from 

loops (hulls) that pierce the piston and the piston head, hut not the cylinder. (11) 

The cylindrical symmetry of the domains in Fig. 1 allows one to rotate q to the positive 
X-axis and perform one angular integral trivially, resulting in 

2 nco noo nco 

£int(cyl.piston) = - — ( / pdp / dz d\^CkQ[{p,^, z) + S)[li{Q)] C \'i^k\)i,(o) ■ 

(12) 

The remaining integrals have support in the region defined by, 

a) 3 i G {1, . . . , f }; z + 2;j > Aa , and 

h) V i G {1, . . . , t;}, {xi + pf + yl < (Ar)^ , and 

c) 3ie{l,...,v}; z + Zi<0 A (x, + p)^ + yf + {z + Zi ~ X^W^f > {XRf . 

(13) 
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Figure 2: On the left is a typical triangulated hull of a Brownian bridge defined by 10^ points. It 
has 286 faces and 145 vertices and the shape of an irregularly cut rhinestone with large triangular 
facettes and intricate corners. The average number of vertices (v) of the convex hull of a Brownian 
bridge given by n points can be read off the left axis of the graph using the upper (red) curve. 
Data points are from numerical simulations. The inset is the formula for the (red) trendline that 
best describes the data. The right axis and lower (green) line give the average CPU-time used to 
compute the hulls on a laptop with a 2GHz processor. Note that the left axis is linear while the 
right axis is logarithmic. 

Condition [T3b) implies that the convex hull pierces the piston at height a, [T3b) that it does 
not pierce the cylinder of radius r and[T3b) that it pierces the cylinder head of radius R> r. 

We first determine the maximum (zmax) and minimum (-Zmin) values of the set {zf, i = 
1, . . . ,f} for a convex hull. Given p G [0, (-Zmax — -Zmin)^/^] and the corresponding minimal 
value of A implied by [T3b). the integration regions of A and z are then found by solving the 
non-linear optimization problem of Eq. f|T3|) using a plane sweep algorithm[3T] that deter- 
mines intersections and endpoints of the quadratic curves in [T3b). Note that the quadratic 
functions of fT3b) decrease monotonically for z < —Zi. Once the piecewise parabolic bound- 
ary of the integration region has been obtained, the z- and A-integrals of Eq.l ITOl) over this 
region are performed analytically. 

The remaining one-dimensional integral over p has to be evaluated numerically. We com- 
pared two methods: i) averaging the results of numerically integrating p (using an adaptive 
algorithm) for every hull separately and ii) averaging (scaled) values of the p-integrand of the 
hulls at gaussian quadrature points and performing the numerical p-integration once only. 
For comparable accuracy both methods require about the same computing time since most 
effort is spent in determining the integration region, rather than in performing the integral: 
for comparable accuracy, the adaptive algorithm on average requires less evaluations of the 
integrand than gaussian integration. 
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To verify the accuracy of the algorithm and estimate systematic errors, we compare (see 
Fig. 3) the numerical results for R ^ r ^ d with the known analytic Casimir energy due 
to a scalar satisfying Dirichlet boundary conditions on two parallel circular plates of area 
S = vrr^ (half the electromagnetic Casimir energy [Ij of this configuration). For R > r ^ a 
the numerical results are compared with the analytical estimates of the asymptotic behavior 
discussed below. 

Piston height a/r 

0.00 0.02 0.04 0.06 0.08 0.10 0.12 




R=r 
R=1.02r 
R=1.05r 
R=1.02r 
(PFA) 

R=oo 



Figure 3: The dimensionless Casimir interaction energy r£int{a) / {he) for the cylindrical cavities of 
radius r with caps of radius i? > r of Fig. 1 as a function of the rescaled height a/r of the moveable 
piston. Dots indicate numerical results (dot size does not represent errors, which are too small to 
show on this logarithmic plot). Solid lines are piston-based Proximity Force Approximation (PFA) 
estimates (note that the PFA overestimates the interaction energy by an order of magnitude for 
R = 1.02r). Dashed curves correspond to the force given by the asymptotic semiclassical estimate 
of Eq. (j2ip . Note the many orders of magnitude between a piston with a flat (bottom, purple) and 
a hemispherical (top, red) cylinder head. 



4 Estimate of the Casimir energy of a cylindrical pis- 
ton with spherical cylinder head for a/r <^ r/R < 1 

4.1 Direct heuristic approach 

The asymptotic contribution to the Casimir energy of the cylindrical piston in Fig. 1 for 
small a/r <^ 1 is dominated by short loops of length 0{a) that satisfy Eq.( fT3i) . Since a 
unit loop generally has extension in 2;-direction of Az = 0{1), the main contribution when 
a/r -C r/R comes from the behavior of the integrand of Eq. lfTOj) at large A ~ Az/a (small 
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f3) and for r/R ^ 1 is concentrated near p ~ Ar ^ Xi,yi for all vertices. In this regime the 
boundary of the integration region in Eq. (fT3|) approaches, 




■mm 



- Z, 



■mm 



p < Xr 



(14) 



where z^^x and z^^i^ are the maximal/minimal 2;- values of the vertices of the unit loop. 
Eq. (HM apparently is invariant under rotations of a hull about the z-axis, but only because 
it was implicitly assumed that for any given A the hull is oriented so that conditions ( fT3l) 
are satisfied for the largest possible value of p. As can be seen from Fig. 4, this implies 
choosing a particular hull vertex as the origin of the hull and orienting the hull optimally 
with respect to the casing by rotating it about its z-axis. Selecting just one particular 
representative of a class of rotated hulls turns out to be a reasonable approximation for 
obtaining the asymptotic behavior when r ~ /2. [When r = R, the fact that the smallest 
possible hull of a given shape must not pierce the cylinder wall while piercing the piston and 
hemispherical cap evidently greatly constrains the orientation of the hull with the greatest 
weight - slightly rotating away from the optimal orientation requires smaller p and A or may 
even make it impossible to satisfy ( fTTl) ]. Note that conditioning on hulls with a particular 
orientation implies that the volume Vso{2) = 27r of rotations about the z-axis should be 
divided out of the rotation invariant probability measure for the hulls. Requirements [T^) 



Figure 4: An optimally oriented convex hull of a unit loop. The hull is rotated about a vertical 
axis to maximize its contribution to i?int(fl) in Eq. (jl2p for small piston height a. Almost all hulls 
have the shape of an irregularly cut rhinestone (the one depicted here has 55 vertices and is that of 
a Brownian bridge defined by 10^ points). There generally is just one optimal relative orientation 
of the hull with the greatest weight. Rotation symmetry about the cylinder axis has already been 
taken into account in Eq. (jl2p . Note that for small a and a hemispherical piston, a less optimal 
orientation leads to a great reduction of the available phase space, since the hull must satisfy 
condition pTj) and pierce the piston and cap but not the cylinder. 
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and [T^ ) imply that Xa — z^^s_^ < z < Xy/R'^ — — a/ (Ai?)^ — — 2;min. The upper bound 



of the p-integration is Xr and its lower bound is y {XRy — (AV-R^ — + A2; — Aa)^ when 

A-R > XV — + — Aa, and zero otherwise. Neglecting terms that manifestly are finite 
in the limit a/r 0, the asymptotic behavior of Si^t for a/r <^ r/R (upon substituting 
= Xa, X = Az/'j) is given by, 



An Vso{2) J„ y_ , 



Ji2-(V-R2_r2+T._a)2 



967r2a 



^ + 1 + ) ((Az)^).,(o) . (15) 



The expectation {{AzY) i-^^q) is related to the Casimir interaction energy of two parallel 
discs of radius r, i.e. the limit r/R<^a/r<^loi the piston in Fig. 1. The integration 
region (fT3l) in this case is readily found and one obtains, 



^int(0~^<f <1) ~ -^^l iAz-Xa)dp'')e,io) 

= -i^(l + ^^^/^))^(^^)')M0)- (16) 

Note that we have not divided by ^50(2) = 27r in Eq. (fT6l) because the contribution from loops 
near the cylinder surface is negligible in this limit and no particular orientation generally is 
much preferred. Comparing to the known Casimir interaction energy due to a scalar field 
for two parallel discs of radius r one has that, 

((A^)')4(o) = ^ , (17) 

a result for unit Brownian bridges that can also be found directly. For a/r -C r/R the 
interaction energy of a Casimir piston in Fig. 1 with a cylinder head of radius i? ~ r thus is 
estimated to asymptotically approach. 



^int(f « i < 1) = (^^^^^ + 1 + 0{a/r) ) . (18) 



Including only statistical errors, the asymptotic numerical data for r < < 1.02r is best 
reproduced by. 



The coefficients of the \/R? — r^/a^ and 1/a terms indeed are comparable and the leading 
0.00326(4)/a2-behavior of the attractive force on the piston with a hemispherical head differs 
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by less than 5% from the estimate obtained from Eq. (fT8|) . A semiclassical determination of 
the asymptotic behavior gives a coefficient that is even closer to the one obtained numerically. 
The semiclassical analysis below furthermore shows that the asymptotic behavior of £iat is 
the same in magnitude but of opposite sign for a scalar field satisfying Neumann boundary 
conditions. 




Figure 5: Some closed classical paths that depend on the piston height a and become arbitrary 
short as a — > 0. a) closed paths with just one reflection, b) paths with two reflections, c) paths with 
three reflections; similar closed paths with more reflections are not shown. As explained in the main 
text only paths of type a) with just one reflection contribute to the asymptotic interaction energy. 
Contributions due to paths b) and c) vanish for Neumann and Dirichlet boundary conditions (in 3 
dimensions). 



4.2 Semiclassical approach 

Semiclassically a Casimir energy diverges for some values of the geometric parameters due to 
closed classical paths that become arbitrary short in this limit. The contribution of periodic 
orbits to £int{a = 0) of a hemispherical piston for a scalar field satisfying Dirichlet boundary 
conditions was obtained in [32]. 

^npr / N ^CTT / TT^ 30-\/2 cosfmvr/ri) \ ^ ,^ , 

C(« = 0) = T^ 1 + ^ + 5^ ; . y : ^ 0.044239... nc/i?. (20) 

\ m=l n=2m+l ^ / 

The length of any periodic orbit of the hemisphere exceeds 2r and the interaction Casimir 
energy due to periodic orbits therefore is finite even at a = 0. However, for Dirichlet (or 
Neumann) boundary conditions there are additional contributions to £^int from non-periodic 
closed classical paths that depend on the position of the piston and become arbitrary short 
as a — i> 0. The closed paths shown in Fig. 5 potentially give a divergent contribution 
to the Casimir energy in the a — limit. The trajectories of Fig. 5b) and Fig.5c) have 
an odd number of conjugate points and do not contribute [28]. Independent of this phase 
cancelation, the length of the trajectories in Fig.5b) and Fig. 5c) increases rapidly with R and 
they therefore do not contribute to the Casimir energy for R — > oo. The paths of Fig. 5b) 
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and 5c) thus cannot possibly explain the heuristic and numerical result that the Casimir 
energy of a hemispherical piston does not vanish for r = i? ^ oo. However, contrary to 
the electromagnetic casejl], the classical trajectory of Fig. 5a) that reflects just once off the 
piston for Dirichlet (D) or Neumann (N) boundary conditions gives rise to the asymptotic 
force, 

he r pdp 
^10^ Jo (a + ^/W^ - - r2)"4 

Only terms that diverge for a — have been retained in this estimate, because other non- 
singular contributions to the force of (^(l/r^) have also not been included. Apart from a 
3% smaller overall strength (fy^^j = 1.033. . .), the force in Eq. fl2Tl) is consistent with the 
asymptotic behavior in Eq. (fT8|) that was obtained heuristically. The semiclassical calculation 
in particular also shows that the asymptotic 1/a^ behavior of the force vanishes like \/B? — r"^ 
a.s R ^ r. For a hemispherical Casimir piston, the force for a ^ r is proportional to 
and (somewhat surprisingly) does not depend on the radius r = R oi casing and cylinder 
head. Note that the semiclassical coefficient l/(967r) ~ 0.00331 ... of this asymptotic force 
is only 2% larger than the coefficient 0.00326(4) found numerically. The 2% discrepancy 
probably is as good an estimate of the systematic errors of the numerical calculation as any 
one can give at this level of accuracy (the statistical error of the calculations is about 1.2%). 
Fig. 6 shows the difference between the numerical results and the asymptotic semiclassical 
energy for the hemispherical Casimir piston in a linear plot. Note that the residual force is 
small and repulsive. 



Fig. 3a) 
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5 Discussion and Conclusions 

The Casimir force on a piston by a massless scalar field satisfying Dirichlet boundary condi- 
tions depends qualitatively on the shape of the casing in the systems of Fig. 1. Whereas the 
attractive force for small piston height a/r <^ 1 is proportional to r^/a'^ for a flat cylinder 
head, it is proportional to \/R? — r'^ /a^ for heads with curvature radius r < R <^ r'^/a and 
for a hemispherical piston with {R = r) is proportional to 1/a^. Note that the asymptotic 
behavior of Eq. fl2T]) implies that one percent deviation in radius from a hemispherical cylin- 
der head doubles the Casimir force on the piston at an elevation a/r ~ 0.1. The fact that 
the attractive Casimir force is greatly reduced for a hemispherical cap is well born out by 
the numerical calculations. Fig. 3 indicates that the Casimir force on the piston depends 
qualitatively on R/r and at small piston height a/r may differ by several orders in magni- 
tude. [Because a larger cap radius implies a greater "average" separation between the piston 
and the cap, the overall change in magnitude is somewhat deceptive, but qualitative changes 
in the Casimir force in this logarithmic plot appear as a change in slope and the reduction 
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Figure 6: The difference AS (a) = Smtia) + hc/{9G7ra) of the numerical- and asymptotic 
semiclassical- estimate for the scalar Casimir interaction energy of a hemispherical piston. The 
energy difference is plotted in dimensionless units against the piston height a in units of the cylin- 
der radius r. Only statistical errors of the numerical calculation are shown. The numerical results 
are those of Fig. 3 for R = r and were obtained from the convex hulls of 10^ unit loops, each given 
by 10^ points. The (red) dot at a/r = is the semiclassical contribution to this difference due to 
periodic orbits given by Eq. (f20l) . 

of the Casimir force in a hemispherical piston clearly is not simply due to the increase in 
average piston height.] 

The interaction Casimir energy defined in Eq.l fTOl) is a priori finite (and does not require 
regularization) when the conditions of Eq.(l5]) are satisfied. Only a portion of the vacuum 
energy is computed that includes all its dependence on the piston height a and determines 
the force on the piston. To relate this interaction vacuum energy of a massless scalar to the 
spectral function of the Laplace operator, leading terms in the high-temperature expansion 
have to be canceled. Here this was achieved by subtracting spectral functions of domains with 
the global characteristics represented by the first few (five in three dimensions) heat kernel 
coefficients. The subtraction procedure in particular implies that finite Casimir energies are 
differences in the vacuum energy of domains with the same topology, same volume, same 
average curvature etc... 

The linear combination of vacuum energies for Si^t in the present examples is shown in 
Fig. 1 - for L ^ r only Sq) depends on the height a of the piston. The Casimir energy due 
to a massless Dirichlet scalar in this case is given by the positive probability measure for 
Brownian paths that satisfy Eq. ffTTl) . This implies an attractive Casimir force on the piston 
for any separation a and cap radius R > r. The same argument implies that the Casimir 
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force by a Dirichlet scalar is attractive for a large class of piston geometries for which a 
condition on the Brownian bridges like Eq. (fTT|) holds - this includes pistons and cylinders 
of arbitrary cross-section and form with caps that are entirely within the cylinder. Cylinder 
heads that extend beyond the cylinder have also been examined pl3]. 




For convex Euclidean domains, the spectral function is related to a probability measure 
on convex surfaces. This measure is indirectly induced by that on Brownian bridges the 
convex surfaces are the hulls of. Convex surfaces in this investigation were constructed from 
the Brownian loops and were not generated directly. It would be very advantageous for the 
accuracy and speed of numerical simulations to directly construct convex surfaces with the 
appropriate measure. 

A semiclassical analysis of the geometries presented here reproduces the numerical results 
for the force on the piston to better than 2% when a/r < 0.1, it in particular gives the 
correct asymptotic behavior of the force for a/r ^ whereas the PFA is inapplicable and 
overestimates the force by an order of magnitude. For a massless scalar satisfying Dirichlet 
boundary conditions, this asymptotic behavior is semiclassically given by the contribution 
of the single closed classical path of Fig. 3a) that reflects just once off the piston. The 
contribution to the Casimir energy of periodic orbits calculated in [32] is finite at a = and 
negligible for a/r <^ r/R. The non-periodic path of Fig. 3a) determines the asymptotic 
behavior of the Casimir energy because its length is 0{a) near the edge of the piston, 
whereas all periodic paths have a length exceeding 2{R — ^R? — r^). For a hemispherical 
cap the attractive Casimir force on the piston —hc/{96'n'a'^) does not depend on the radius 
r = R of the cylinder and cap and mimics the electrostatic force on a metallic piston a 
charge /{he) = l/(487r) ~ 1/150.8 < a^M at the center of the spherical cap would exert. 
Although the same as the force due to a charge, the Casimir force here is inherently non- 
local with an energy density that is greatest near the intersection of cylinder and cap. The 
(non-local) nature of this force becomes drastically apparent if one lets the cylinder- and 
cap- radius tend to infinity while keeping the height a of the piston constant: the Casimir 
force on the piston (locally a flat plate) remains unchanged in this limit and the direction 
of the force would give information about the position of the piston relative to the cap at 
infinity. This paradoxical situation arises only for the Casimir energy of a massless scalar 
field satisfying Dirichlet (Neumann) boundary conditions. Semiclassically the asymptotic 
force due to a scalar field satisfying Neumann boundary conditions in Eq. lpT!) is of the same 
magnitude but of opposite sign. [In the method of images the "mirror source" for Neumann 
boundary conditions at the piston is of the same sign as the source, whereas it is of opposite 
sign for Dirichlet boundary conditions]. 



The semiclassical calculation in[32] suggests that the electromagnetic Casimir force on a 

piston with hemispherical cap vanishes for r —>■ oo and is weakly repulsive. Because closed 
paths with an odd number of reflections vanish in the electromagnetic case[l], one does not 
expect the electromagnetic interaction energy of a hemispherical piston to diverge [M] for 
a — >■ 0. Indeed, if the electromagnetic Casimir force for small a/r is (9(l/r^), one cannot 
deduce the relative orientation and position of cap and piston for r —>■ oo. [One does not 
encounter this problem with locality for two fiat plates: for a <^ r the distance between 
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the plates and their relative orientation can always be found locally and does not require 
knowledge of the situation at distances 0{r) from the center.] In contrast to Dirichlet 
boundary conditions, a repulsive force on a piston in a metallic cylinder does not violate 
reflection positivity, because the metallic cylinder induces correlations in the electromagnetic 
fluctuations on either side of the piston. That vacuum forces in some geometries are small 
or at least comparable to electrostatic forces down to very small separations could have 
implications for precision measurements of forces and the design of Micro-Electro- Mechanical 
devices. 
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